Interference of Fock states in a single measurement 
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We study analytically the structure of an arbitrary order correlation function for a pair of Fock 
states and prove without any approximations that in a single measurement of particle positions 
interference effects must occur as experimentally observed with Bose-Einstein condensates. We also 
show that the noise level present in the statistics is slightly lower than for a respective measurement 
of phase states. 
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Paul Dirac in his famous textbook on quantum me- 
chanics [H describes photon interference in the following 
way. Suppose we have a beam of light consisting of a large 
number of photons split up into two components of equal 
intensity (■■■)■ If the two components are now made to 
interfere, we should require a photon in one component 
to be able to interfere with one in the other. Sometimes 
these two photons would have to annihilate one another 
and other times they would have to produce four pho- 
tons. This would contradict the conservation of energy 
(...). Each photon then interferes only with itself. Inter- 
ference between two different photons never occurs. This 
view has been criticized |2j for applying only to the states 
of a definite number of particles. For example two inde- 
pendent sources of coherent states which have no defi- 
nite energy can generate an interference pattern without 
question. 

Dirac's argument, however, may seem to apply at least 
to the particle-number Fock states that do have a definite 
energy value. Another reason why Fock states seem inca- 
pable of interfering is that they do not have well defined 
relative phase. And the last reason - a direct calcula- 
tion of the first-order correlation function for two Fock 
states does not reveal any interference properties. Unfor- 
tunately, these three very attractive arguments fail. 

The first beautiful experimental example of the Fock 
states interference has been accomplished with Bose- 
Einstein condensates that can be thought of as particle 
number states. In the double-condensate experiment by 
Andrews et al. Q the authors proved the existence of 
interference fringes in the measurement of positions of 
condensate atoms. Results of a similar experiment have 
been recently reported in Ref. 0|. How is it possible? 

The reason is that the first-order correlation function is 
attributed to an average (over many realizations) density 
of particles, while in the experiments 0, 0] we deal with 
the results of a single measurement. Quantum mechanics 
cannot predict the exact result of a single measurement 
- it only predicts average values of certain observables 
or a probability of a definite result of a single experi- 
ment. But since we are dealing with a huge number of 
particles, how about using this many particle probability 
distribution to predict a typical particle density profile in 
a single experiment [H|? Apparently considering only the 



first-order correlation function is not enough to guess the 
typical density shape and one needs to take into account 
also higher-order correlation functions. 

This issue has been first addressed in a beautiful work 
of Javanainen and Yoo [6] where the authors apply nu- 
merical analysis to the studies of the structure of many 
particle probability density distributions. In this and the 
subsequent numerical experiment Q exploiting the laws 
of quantum mechanics the authors show that two Fock 
states can indeed reveal an interference pattern in a single 
measurement. Obviously after averaging out over many 
realizations of the numerical experiment the interference 
effects disappear as expected. 

In this Letter we analytically study a nature of the 
high-order correlation functions to show directly from 
their mathematical structure the existence of interference 
effects in a single interference measurement of two Fock 
states. In our analysis we do not use any approximations, 
as the previous authors who attempted to prove this 
result analytically [3, |9| assuming orthogonality of the 
phase states. This approximation is questionable when 
the number of particles measured is of the order of the 
total number of particles of the system. We also show an 
interesting and unintuitive property of the noise present 
in the interference pattern. The noise level turns out to 
be slightly lower than in the case of multiple drawn po- 
sitions with a probability distribution equal to the inter- 
ference pattern. We have shown in a numerical test that 
the difference is very small and probably out of reach of 
any experimental observation. Our result, however, gives 
an interesting insight into the structure of the high-order 
correlation function. 

Consider a set of d identical ideal detectors capable of 
counting particles. Let the surface L of the z-th detec- 
tor placed at the position cc, be described by the char- 
acteristic function x( x ~ x i) an d the annihilation opera- 
tor Ai associated with the mode L~^xi x — x i)- Let us 
assume, that the detectors are spatially separated, i.e. 
J x( x — x i)x( x ~ x j)dx = L5ij. We calculate an average 
product of the particle counts from all d detectors [Toj | : 



I(xi,...,x d ) = (A\Ai 



■A\A d ) 



(1) 



The set of particles measured by the detectors is de- 
scribed by the field operator ^>(x). We assume, that the 
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occupied modes are slowly-varying in comparison to the 
size L of the detectors: 

Ai = J -jjX{x - Xi)$(x)dx w VZ4f(xi). (2) 

From the above formula follows a connection between 
average product of detector counts with the d-order cor- 
relation function: 

I( Xl , ...,x d ) = L d (¥( Xl ) . . . ¥(x d )*( Xl ) . . . 

(3) 



If we assume that the detectors' size L is so small that 
each of them detects, on average, much less than a single 
particle then the average product of the particle counts 
I(xi, . . . ,x d ) can be identified with a probability of de- 
tection of exactly one particle by each detector. Thus the 
probability density g of localizing the first particle at the 
position xi, the second particle at X2, etc., equals: 



,x d ) 



(N-d)\ 
M 



x (&(x 1 )...¥(x d )*(x 1 )...i>(x d )), 



(4) 



where N is the total number of particles. Let us notice, that the probability density Q is defined only for the states of 
a definite number of particles. This approach allows one to interpret the physical meaning of the correlation function 
of the order d in two ways. On the one hand it is proportional to the average product of particle counts of d detectors, 
on the other hand it is related to the probability density of localizing exactly one particle by each of d very small 
detectors. 

As long as the detectors are spatially separated an ordering of the field operators in the expression ([3]) and Q is 
defined up to the commutation relation [&(x), = S(x — y). If one wants to continuously extend the expressions 

to the case of X4 = Xj for i 7^ j then the field operators must be ordered normally. 

Consider a two-mode quantum state |n, N — n) with the first mode defined by an arbitrary function u(x) and the 
second orthogonal mode by w(x) for x E [0,1]. From the expression |(4|) we calculate the probability distribution of 
localizing all the N particles at positions x%, x%, . . . , xjf. 



Q\n,N-n){ x l, ■ ■ ■ 1 Xn) 



Z2 U ( X V(.1)) ■ ■ ■ u {x v{n) )w{x v(n+1) ) . ..w{x v(N) ) 



(5) 



where we sum up over permutations V of an TV-element set excluding the non-trivial permutations acting separately 
on the first n elements of the set and the last N — n elements. We will consider the case N — 2n, when exactly n 
particles occupy each mode. Using the formula |(4]) we find that the probability density of detecting d of 2n particles 
at the positions X\, x%, . . . , x d can be expressed with the probability densities for the asymmetric states |(5]): 



Q\n,n)(x\ 



,X d ) = 



-1 d 

E 



0(n- j)Q(n-d + j) 



2n-d 
n — d + j 



Q\j,d-j)(xi, ■ ■ -,x d ), 



(6) 



where <3(x) is the Heaviside's theta function. Binomial coefficients 



showing up in the above expression 



Kk) ~ k\(n-k)\ 

are for n^$> 1 bell-shaped functions of k centered around = f and with a dispersion equal to . We see that the 
coefficients ( „ 2n 7 + ■) and ( .) attain their maxima for the same value j = | , but they are characterized by the different 



\n-d+jj \j 
dispersions of the variable j: 



and respectively. 
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Consider a special case of the probability density ([6]), with only a small fraction of all particles being measured, 
d <a In this case the distribution (n™ d +j) is much wider than (*) and we can replace the former with its maximum 
value. In this case also the Heaviside's thetas are equal to the unity and we can skip them. As a result the expression 
((6]) can be written in the following form: 



Q\n,n){x\ 



,X d ) 



7=1 V7/ J ~ 7r 1=1 



(7) 



We have managed to express the low-order correlation function for the highly occupied state \n, n) in an ele- 
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gant form of an integral over some positive expression. 



A similar result has been shown in Refs. [8|, [y], how- 
ever the authors omit the fact that they actually prove 
it only for d <C n because of the limited validity of the 
approximations used. These approximations are highly 
questionable when the number of particles measured is 
of the order of the total number of particles d ~ 2n, 
therefore we are going to prove all the properties of the 
high-order correlation functions with no approximations 
whatsoever. 

It turns out, that the Eq. J7]) tells a lot about a re- 
sult of a single measurement of positions of d particles. 
According to Born's probabilistic interpretation of quan- 
tum mechanics a result of such measurement - the set 
of measured positions xi,...,Xd, corresponds to a re- 
sult of a single drawing with the probability density 
Q\n t n)(xi, ■ ■ ■ ,Xd)- Let us try to predict the result of 
such drawing using the following simple lemma based on 
Bayes' theorem: 

Lemma. If N -dimensional probability density g 
can be represented in the form q{x\,...,xn) = 
J d£p(£)g(xx, . . . , xn\£), where p is a one- dimensional 
probability density and q is an N -dimensional conditional 
probability distribution (likelihood) then drawing a set of 
random variables (x±, . . . , xn) with the probability g is 
equivalent to drawing a random variable £ with the den- 
sity p, and then drawing the set of random variables 
(x\, . . . , xn) with the density q for the chosen £. 

Proof. Equivalence of both densities can be shown by 
proving the equality of arbitrary moments of the distribu- 
tions. We will use an elementary theorem about changing 
the order of integrals. An arbitrary moment for the sec- 
ond distribution reads: 



y"d£p(£) J dxi...dx N x* 1 



' N 



fq(xi,...,x N \£) 



dxi... dx N xl 1 ... x k N N I d£p(£)q(x 



and it is equal to the same moment for the distribution 
g. As we know, the values of all the moments uniquely 
determine the probability distribution. QED. 

It follows that the result of a single draw with the prob- 
ability density (7]) can be achieved by a preliminary draw 
of the parameter <f> with a flat distribution, and then 
by drawing positions of particles according to the sep- 
arable density n»=i 2 I^C^t) + e^w(xi)\ . The second 
draw yields positions centered around maxima of the one- 
dimensional function \ \u{x) + e l ^w(x)\ . If we assume 
u(x) = w*(x) = e t7rx , x G [0,1] then every single mea- 
surement reveals the interference fringes with maxima 
located randomly each time somewhere else. The mean- 
ing of the interference fringes can be made more precise 
in the following way. Suppose that the whole space of 
possible particle positions is divided into D small areas 
of equal length and we examine how many of the first d 



particles enter each of these areas in a single measure- 
ment [ll|]. Each area is tightly covered by a set of small 



detectors constituting, so to say, a single super-detector. 
We look at the histograms of the count statistics of the 
single measurement - if the sizes of the considered areas 
are such that each of them swallows on average a large 
number of particles, then each histogram should repro- 
duce the function \ \ u(x) + e l( ^w(x)\ for some cf>. 

Let us notice that the last expression in the formula 
([7]) defines a hidden variables model with the role of hid- 
den parameter played by <fi. Therefore the uncertainty 
of the phase 4> attributed to the single measurement of 
the small portion d of all particles must be of a classi- 
cal nature. Although one can establish a link between 
a spin-| formalism and parameter <\> our observation in- 
dicates that the considered type of measurement cannot 
lead to violation of Bell's inequalities. 

We have just shown that for the single measurement 
of the relatively small number d of particles belonging 
to the state \n, n) one observes the interference fringes. 
Therefore it is natural to ask about the result of a simi- 
lar measurement of all the 2n particles. Below we show 
that the larger number of particles is being measured, the 
fringes of even higher quality are observed. This agrees 
with the numerical test [3] and the methodology of the 
experiment fill ]. 

The proof is the following. According to our lemma 
drawing d positions described by the probability distri- 
bution (7]) can be achieved by drawing parameter </>, and 
then drawing positions with the conditional probabil- 
ity distribution n<=i 3 \ u i x i) + e % ^w(xi)\ . However, ac- 
cording to the expression ([7]) another equivalent method 
of drawing exists and is based on drawing first the pa- 
rameter j described by the distribution 2~ d (^) and then 
drawing positions of the particles with the probability 
density Q\j d-j) i x ii • • ■ i x <i) given by the analytic formula 
([!]). Obviously, both methods of drawing lead to the 
same result which, as we know, reveals the interference 
patterns of the known shape. The second equality in J7]) 
indicates that the results of drawing of positions with 
the probability distribution Qu^-j) ( x ii ■ ■ ■ > x d) f° r the 
parameter j differing from | by not more than a few dis- 
persion lengths ^ must reveal the interference effects 
every time. Independently of the method of drawing 
each random histogram will vary from the ideal shape 
\ \u(x) + e l</3 w(x)\ because of statistical fluctuations. 
Let us introduce the following measure of these fluctu- 
ations defined for an arbitrary result of the single draw- 
ing. Let the number of counts of the z-th super-detector 
placed at Xi be denoted with n,-. For the histogram of 
results {n>i} we define the following quantity: 



X 



inf 



D 



i=l 



d 

2D 



(8) 



where D is the number of the super-detectors. The 
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above expression averaged out over many realizations 
X 2 we will call noise. This noise depends only on the 
number of super-detectors and the probability distribu- 
tion g, or equivalently on the quantum state g and the 
parameters d and D, which we denote as x 2 {Qi d, D). 
Therefore the better the histograms reproduce the shape 
\ \u{x) + e l ^w{x) | (for some <f>) the lower the value of 
noise x 2 ■ From the last equality in the formula (J7j) we 
get: 

]T2- d ( d V(|j, d-j),d,D) = ^(\d}^d,D), (9) 
j=i V- 7 / 

where \d)$ is so called phase state of d particles occupy- 
ing the same mode A= [u(x) + e l ^w{x)\ . We have used 

the fact that the quantity X 2 (\d) d, D) cannot depend 
on the selection of <j> and it determines the level of noise 
for the histogram of particle positions drawed one by one 
with the probability density \ \u(x) + e l ^w(x)\ . Equa- 
tion (J9j) indicates that the noise level x 2 {\d)<f>, d, D) is 
equal to an average noise level for the states \j, d — j) 
with the weights equal to 2~ d (^). 

Unfortunately the level of noise averaged out over all 
states \ j, d — j) does not uniquely determine the value of 
noise X 2 (|ij d — j), d, D) for the particular j. However we 
can use a natural assumption that the interference effects 
disappear for the asymmetric states \j, d — j). In the ex- 
treme but highly improbable example of the state \d, 0) or 
|0, d) the interference will be obviously completely absent. 
To be more specific, we assume that X 2 (|i, d — j), n, D) is 
a monotonically increasing function of \j — ||, Accord- 
ing to this assumption and the equation |(9|) we anticipate 
the following inequality to hold: 

5?(|d/2, d/2), d, D) < x*(\d)t,d, D), (10) 

which completes the proof. 

We have investigated validity of this inequality by com- 
paring the noise in numerically drawn histograms for 
the state |d/2,d/2) and \d)<f, but the observed differ- 
ence did not exceed the level of statistical error. This 
means that the difference between x 2 (|d/2, d/2), d, D) 
and x 2 (\d)<p,d,D) is very small, which reflects the fact 
that the probability of drawing the highly asymmetric 
state in (7]) is negligible. The non-intuitive inequality 
(fTO)) . although very weak, must be an interesting signa- 
ture of non-trivial spatial correlations present within the 
mathematical structure of the Fock states. 

The inequality (fTUf can be seen also from the struc- 
ture of the analytic expression ((6]). Let us notice that 
when the number d of the drawed particles approaches 
its maximum value 2n then the width of the distribu- 
tion ( n 2 "rf_i) equal to V2 %~ d rapidly shrinks. It follows 
that the more particles we measure, the more symmetric 
states (which more likely contribute to the interference) 
are being chosen for the drawing of the positions. 



Our last conclusion is that in the limit of d ^> 1, the 
probability distribution 2~ d (^) from the Eq. ([7]) becomes 

relatively narrow as ^ <C d and only the states \j, d — j) 
that are almost symmetric will be chosen for the drawing 
of the particle positions. Therefore in the large particle 
number limit all quantities that weakly depend on the 
asymmetry of the state will reproduce the results ob- 
tained for the phase states \d)j,. 

We have proven the existence of the interference effects 
by studying the structure of the high-order correlation 
functions for the Fock states. It is also clear that these 
effects will disappear after averaging out over many rep- 
etitions of the measurement. This result is, however, an 
immediate consequence of the Bogoliubov method which 
assumes ad hoc that one can replace the field operator of a 
single condensate by a classical wave with small quantum 
corrections: * « vNe 1 ^ + 5$ and neglecting the latter. 
Our analysis allows one to attribute the arbitrarily cho- 
sen phase 4> in the Bogoliubov method with the param- 
eter <p from the equation ([7]) spontaneously induced in a 
single measurement. In this interpretation breaking the 
phase-space symmetry of the Fock states by using the 
Bogoliubov method corresponds to replacing the strict 
expressions given by Eq. J6]) with their approximations 
©. 
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